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ABSTRACT 

O 

^ Extrasolar terrestrial planets with the potential to host life might have large 

obliquities or be subject to strong obliquity variations. We revisit the habit ability 
of oblique planets with an energy balance climate model (EBM) allowing for 
Q dynamical transitions to ice-covered snowball states as a result of ice-albedo 

^ feedback. Despite the great simplicity of our EBM, it captures reasonably well 

the seasonal cycle of global energetic fluxes at Earth's surface. It also performs 
C — satisfactorily against a full-physics climate model of a highly oblique Earth-like 

^ planet, in an unusual regime of circulation dominated by heat transport from the 

O poles to the equator. Climates on oblique terrestrial planets can violate global 

^ ^ radiative balance through much of their seasonal cycle, which limits the usefulness 

^ of simple radiative equilibrium arguments. High obliquity planets have severe 

climates, with large amplitude seasonal variations, but they are not necessarily 
more prone to global snowball transitions than low obliquity planets. We find 
that terrestrial planets with massive CO2 atmospheres, typically expected in the 
outer regions of habitable zones, can also be subject to such dynamical snowball 
transitions. Some of the snowball climates investigated for C02-rich atmospheres 
experience partial atmospheric collapse. Since long-term CO2 atmospheric build- 
up acts as a climatic thermostat for habitable planets, partial CO2 collapse could 
limit the habitability of such planets. A terrestrial planet's habitability may thus 
depend sensitively on its short-term climatic stability. 



Subject headings: astrobiology - planetary systems - radiative transfer 
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1. Introduction 



The Earth's obhquity is remarkably stable: the angle between the spin-axis and the normal 
to the orbital plane varies by no more than a few degrees from its present value of ^ 23.5°. 
This stability is maintained by torque from the Moon ( Laskar et al.||1993 ; Neron de Surgy & 
Laskar 1997). Even within our own Solar System, though, the obliquity of other terrestrial 



planets has varied significantly more; the analysis of Laskar & Robutel (1993) indicates that 
Mars' obliquity exhibits chaotic variations between ^ 0° and ^ 60°. 



How does climate depend on obliquity and its possible variations in time? How does the 
range of orbital radii around a star at which a planet could support water-based life depend 
on the planet's obliquity? Has the stability of Earth's obliquity made it a more climatically 
hospitable home? The answers to these questions will be important to evaluate the fraction 
of stars that have potentially habitable planets. There are now more than 300 extrasolar 
planets knownQ several of which are close to the terrestrial regime with masses less than 10 
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Mayor et al. 2008). COROT^ which has already launched, and Kepler^ scheduled to launch 
in less than one year, are dedicated space-based transit-detecting observatories that will 
monitor a large number of stars to detect the small decreases in stellar fiux that occur when 



terrestrial planets cross in front of their host stars (Baglin 2003; Borucki et al. 2003, 2007). 
These missions are expected to multiply by perhaps several hundredfold or more the number 
of known terrestrial planets, depending on the distribution of such planets around solar- 



type stars (Borucki et al. 2007, 2003; Basri et al. 2005; although see revised predictions in 



Beatty & Gaudi 2008). NASA and ESA have plans for ambitious future missions to obtain 
spectra of nearby Earth-like planets in the hope that they would reveal the first unambiguous 
signatures of life on a remote world: NASA's Terrestrial Planet Finder and ESA's Darwin 



(Leger & Herbst 2007). The design of such observatories, and the urgency with which they 
will be built and deployed, will depend on the habit ability potential of terrestrial planets 
that will be found in the next 5-10 years. 



Over the last 50 years, various authors have addressed how to predict the way in which 



terrestrial planet habitability depends on star-planet distance (see Kasting & Catling 2003 



for a recent review) . Several of the important initial calculations predated the first discoveries 



of extrasolar planets, including Dole (1964), Hart (1979), and the seminal work of Kasting 



et al.| ( |1993D . [Selsis et al.| (|2007|), |von Bloh et al.| ( |2008[ ), and [Barnes et al.| ([2008|) have 



reconsidered habitability in light of recent exoplanetary detections. Williams et al. (1996), 



Wilhams k Kasting] ( |1997| hereafter WK97) and [Williams fc Pollard] ( [20031 ) have tackled 
precisely the questions relating to obliquity posed above, and have concluded that variations 



^See http: //exoplanet.eu/] 
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in obliquity do not necessarily render a planet non-habitable (see also Hunt 1982, Williams 
T988a||bl IQglesby fc Qgg|[T998l [Chandler fc "Sohl|[2000l and | Jenkins] [20001 [200T| [20031 in the 
context of Earth's paleoclimate studies). 



Here we seek to generalize these analyses to model planets that are less close analogs to 



Earth than have been considered previously. In Spiegel et al. (2008, hereafter SMS08), we 



examined how regionally and temporally habitable climates are affected by variations in the 
efficiency of latitudinal heat transport on a planet, and by variations in the ocean fraction. 
Importantly, we found that otherwise habitable Earth-like terrestrial planets can be subject 
to dynamical climate transitions into globally- frozen snowball states. Since it is not trivial 



to escape a snowball state (e.g., Pierrehumbert 2005) and such globally-frozen climates may 



have profound influences on the development or existence of life (e.g., Hoffman & Schrag 
2002D, identifying the likelihood of such transitions on terrestrial exoplanets should be central 
to any habit ability assessment. Following in the footsteps of our first analysis, here we focus 
on obliquity and consider the influence on habitability of several planetary attributes a priori 
unknown for exoplanets, such as the efficiency of latitudinal heat transport and the land- 
ocean distribution. 



The remainder of this paper is structured as follows: In §[2] we describe the energy balance 
climate model we use. In § [3] discuss several validation tests in which our model performs 
well enough to give us some confidence in its behavior for conditions that differ from those 
found on Earth. In § [4] we examine the influence on regional and seasonal habitability of 
various excursions from Earth-like conditions. Finally, we conclude in § [5} 



2. Model 



In order to describe the surface temperature and its evolution on a terrestrial planet, we 
use a 1-dimensional time-dependent energy balance model based on a diffusion equation for 
latitudinal heat transport. This type of model has been used in previous investigations of 
habitable climates (WK97, SMS08), in modeling Martian climate under changes in forcing 
dNakamura fc Tajika|[2002l [2003[ ), and in studies of the Earth's climate ( [North et aL]|1981 



and references therein). 

Our model is based on the following prognostic equation for the planetary surface tempera- 
ture (as described in SMS08): 

C^-#ffl(l-.^)^U/ = S(l-A). (1) 



dt dx \ dx 

In this equation, x = sin A is the sine of latitude A, T is the temperature, C is the effective 
heat capacity of the surface layer, D is the diffusion coefficient that determines the efficiency 
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of latitudinal redistribution of heat, / is the infrared emission function (energy sink), S is the 
diurnally averaged insolation function (energy source) and A is the albedo. Our formalism 
for insolation on oblique planets follows that of WK97. In the above equation, C, /, and 
A may be functions of T, x, t, and possibly other relevant parameters. 

Our prescriptions for the functions (7, /, and A follow SMS08 and are largely borrowed 
from WK97 or the existing geophysical literature on similar energy balance models (EBMs). 
For simplicity and flexibility, many of our models use very simple, physically motivated 
prescriptions. As described in SMS08, we flnd that an infrared cooling function of the form 

/[r] = arV(i + (3/4)riK[r]) (2) 



(i.e., a one-zone model combined with a simple Eddington transfer approximation; e.g., Shu 



1982), reproduces the greenhouse effect on Earth reasonably well for our purposes. Here, 
TiR represents the opacity of the atmosphere to long wavelength infrared radiation. SMS08 
describes three pairs of infrared radiation functions and albedo functions. In this analysis, we 
will use two of the three: (/2, ^2), which gives the closest match to the Earth's annual average 
temperature distribution, and (/3, ^43), which uses the standard linearized cooling function of 



North & Coakley (1979). The (I^A) functions are detailed in Table [T] The albedo functions 



A2 and As are constant and low (^ 0.3) for high temperatures, constant and high (^ 0.7) 
for low temperatures (to represent the high albedo of snow and ice), and vary smoothly in 
between. One other / function that we use is the one proposed in WK97, derived from full 
radiative-convective calculations, here denoted /wk97- This function includes the detailed 
influence of CO2 (and implicitely H2O) atmospheric content on radiative cooling. For (7, 
we assume various conflgurations of land and ocean, in each conflguration using the same 
land, ocean, and ice partial C values as WK97. Finally, we adopt a diffusion coefficient for 
latitudinal heat transport L^fid = 5.394 x 10^ erg cm~^ s~-^ K~-^ x {Qp/Q^)~'^^ as described 
in SMS08, where Qp is the angular spin frequency of the model planet and is that of 
the Earth. As explained in SMS08, this scaling largely oversimplifies the complexity of 



atmospheric transport expected for planets at different rotation rates (e.g., del Genio et al. 



1993; del Genio & Zhou 1996). 



Equation ([T]) is solved as described in SMS08, on a grid uniformly spaced in latitude (1.25° 
resolution elements, found to be sufficient from convergence tests). We again choose "hot 
start" (T > 350 K) initial conditions, to minimize the likelihood that models will undergo a 
dynamical transition to fully ice-covered (snowball) states from which they cannot recover 
because of ice-albedo feedback. In that sense, our results on snowball states are conservative. 

To summarize, we make the following assumptions in the models presented below: 



1. Heating/Cooling - The heating and cooling functions are given by the diurnally aver- 
aged insolation from a sun-like (1 M©, 1 Lq) star, with albedo and insolation functions 
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described above and in Table [H 

2. Latitudinal Heat Transport - We test the influence on climate of three different ef- 
flciencies of latitudinal heat transport, within the diffusion equation approximation: 
an Earth-like diffusion coefficient, and diffusion coefficients scaled down and up by a 
factor of 9 (which correspond to 8-hour and 72-hour rotation according to the above 
D oc scaling). 

3. Ocean Coverage - We vary both the fraction and the distribution of ocean coverage. 
For ocean fraction, we present a series of models with Earth-like (30%:70%) land:ocean 
fraction, and another series of models that represent a desert-world, with a 90%: 10% 
land:ocean fraction. For ocean distribution, we present models in which there is a 
uniform distribution (in every latitude band) of land and ocean, and others in which 
the land-mass is a single continent centered on the North Pole, while the rest of the 
planet is covered with ocean. 

4. Initial Conditions - As described in SMS08, the models all have a hot-start initial 
condition, with a uniform surface temperature of at least 350 K, to minimize the 
chances of ending up in a globally-frozen snowball state owing solely to the choice of 
initial conditions. Time begins at the Northern winter solstice. 



3. Model Validation 



In SMS08, we verified that our "fiducial" model (70% ocean; /2, A2 cooling- albedo functions) 
at 1 AU predicts temperatures that match the Earth's actual temperature distribution at all 
latitudes that are not significantly affected by Antarctica (i.e., north of 60° S or so). This 
indicates that the model accounts for the overall (annual) planetary energy balance reason- 
ably well. Another obvious test is whether the model correctly predicts the monthly energy 
fiuxes that together go into the overall balance. Because our current investigation tests the 
infiuence of obliquity on climate, and obliquity is the primary driver of the Earth's seasons, 
verifying the seasonal predictions of our model, given Earth-like conditions, is particularly 
relevant. 

The diffusion equation model is a statement of conservation of energy. By definition, after 
vertical integration for a thin atmosphere with dominant surface processes, 

C^ = ^ (3) 



where a is the energy surface density (internal energy per unit surface area on the globe). 
The diffusion equation, therefore, says that the rate of change of internal energy at a given 
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point equals the sources of energy (insolation), minus the sinks (infrared radiation), minus 
whatever energy flows away from the point under consideration. 

Figure [l] presents a comparison between the annually averaged fluxes of incoming and outgo- 
ing radiative energy in the flducial model with the corresponding fluxes on Earth, taken from 



NASA's Earth Radiation Budget Experiment (ERBE) in the mid-1980s (Barkstrom et al. 



1990| )P1 While our model does not capture the full shape of the Earth's cooling and heating 



functions - in particular, the annually averaged model heating function is a bit below the 
Earth's at the poles ^ still, both cooling and heating fluxes are within 10% of the Earth's 
over most of the planet's surface. 

Figure [2] oflFers an even more compelling validation. In this flgure, each of the 12 panels shows 
solar (i.e., heating), terrestrial (i.e., cooling), and net (solar minus terrestrial) radiative fluxes 
as functions of latitude, averaged over one month. Not only are our annually averaged cooling 
and heating functions in reasonable agreement with Earth's, as per Figure [T} but furthermore 
the temporal variability of radiative fluxes in our model is similar to that of Earth. 

For example, at the Northern winter solstice (upper left panel of Fig. [2]), the model heating 
curve closely traces that of the Earth. It peaks at a somewhat more Southern latitude than 
the Earth's does, but is within 10% of the Earth's at all latitudes north of 60° S. As the 
months advance, the concordance between the model heating curve and the Earth's heating 
curve increases, until there is maximum agreement (within 10% at all atitudes) at the equinox 
(''Solstice+03"). Then, by the next solstice, the curves agree to within 10% at all latitudes 
South of roughly 60° N. In a comparison of the cooling curves, the model shows even greater 
agreement with the data. In a majority of months, these two curves are within 10% of each 
other at all latitudes. 



Interestingly, the month-by-month variations in model heating and model cooling lead to 
a net heating curve (heating minus cooling) that predicts some detailed features actually 
seen in the Earth's net heating budget. Notice, for instance, the slight upward turn of the 
net heating curves of both the model and the Earth near the North Pole, at and around 
the Northern winter solstice. A similar feature is seen in both curves (though with slightly 
less impressive detailed agreement) near the South Pole, at and around the Southern winter 
solstice. These comparisons establish that our climate model exhibits reasonable regional 
and seasonal variability of not just temperature but also incoming and outgoing radiative 
energy fluxes. 

Another way to consider seasonal variations of heating and cooling fluxes is to look at 



^The ERBE satellite measured short- wavelength, or incoming, flux as that from 0.2 /im to 4.5 /am. Long- 
wavelength, or outgoing, flux was defined as all other flux within the bolometric range of the instrument. 
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the globe-average of each with respect to time. Figure |3] presents a comparison of these 
fluxes, for our flducial model and the Earth itself. The bottom panel of this flgure shows 
the net heating flux as a function of time of year, measured in fraction of a year from 
the Northern winter solstice. Earth's net heating flux varies by about 5% with respect 
to the cooling flux, while our model's varies somewhat less. The heating function for the 
Earth exceeds the cooling function during Northern winter for two main reasons]^ First, the 
nonzero eccentricity (e ^ 0.0167) of the Earth's orbit places its perihelion - which occurs 
during Northern winter ^ approximately 3.4% closer to the Sun than its aphelion. This is 
responsible for ~ 7% of the annual ~ 10% annual variation in net heating flux. Our flducial 
model, on the other hand, assumes zero eccentricity. Another contributing factor is that the 
Earth's oceans on average absorb somewhat more insolation than the land, and the Southern 
hemisphere - which faces the Sun during Northern winter - has greater ocean coverage than 
the Northern hemisphere. To within 5%, however, the Earth remains in global radiative 
balance throughout its seasonal cycle. 

So far we have considered the radiative fluxes, but what about diffusive energy flux? We 
may combine equation ([T]) with equation ^ to produce: 



cos' A^rj = (sources - sinks)^^^^^^ ^^^^ , (4) 



da d ( ^ o . d 

dt dx 



where we have substituted cos' A for (1 — x'). Comparing equation ([i]) to a diffusion equa- 
tion in spherical coordinates, and accounting for vertical integration, shows that Fx = 
RDco8X{dT/dx) is the rate of latitudinal energy transport per unit longitudinal length. 
The total rate of meridional diffusive heat transport (i.e. energy crossing a given latitude 
circle per unit time) therefore is 

dT 

Tx = 27rRco8XFx = 27ri?'i^ cos' A— . (5) 

ox 



Figure |4] shows proflles of this diffusive heat transport rate in our flducial model, at Earth- 
like 23.5° obliquity and at extreme 90° obliquity. In the Earth-like conflguration, heat flows 
from the equator toward the poles. In the highly oblique conflguration, however, heat flows 
in the other direction, from the poles to the equator (in an annually averaged sense). For 
comparison, Williams & Pollard (2003) present a full general circulation model (GCM) of 
an Earth-like planet at Earth-like and higher obliquity. Figure 2 of that paper shows the 
meridional heat flux within their models for 23.5° obliquity and 85° obliquity, and the results 
are strikingly similar to ours. At 23.5° obliquity, our model's diffusive flux is very close to 
that of the GCM. At high obliquity, the flux in our model remains within ^ 30% of that in 



^Note that the cooling function, which traces surface temperatures, varies less through the seasonal cycle 
than the heating function. 
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the GCM (from visual inspection), at all latitudes. This reasonable concordance indicates 
that the treatment of heat transport within our model, despite being very simple, is still 
likely to remain useful as a representation of heat transport in less-Earth-like conditions. 
We emphasize that it is a nontrivial point that this entirely different regime of transport 
should remain well-captured by a diffusion approximation. 



4. Study of Habitability 

For model planets with 23.5° obliquity on a circular orbit at 1 AU, both pairs of infrared 
cooling functions and albedo functions presented in Table [T] are reasonably good matches 
for the Earth's current climate, as measured by latitudinally averaged temperatures, with 
a somehwat better fit with (/2,yl2)- This gives us some confidence that these functions are 
useful guides as to how the climate might respond under different forcing conditions. In 
this investigation, we consider how variations in intrinsic planetary characteristics combine 
with the changes in insolation and year-length at various orbital radii to map the zone of 
regionally habitable climates on planets with various obliquities. 

We follow SMS08 in saying that, at a given time, a part of a planet is habitable if its surface 
temperature is between 273 K and 373 K, corresponding to the freezing and boiling points of 
pure water at 1 Atm pressure. This criterion may be criticized for several reasons discussed 
in SMS08 and references therein, but it provides a reasonable starting point for making 
numerical investigations. We will frequently quantify habitability of pseudo-Earths with the 
temporal habitability fraction, /timel^^,^], where a is orbital semimajor axis, A is latitude, 
and /time is the fraction of the year that the point in parameter space specified by (a. A) 
spends in the habitable temperature range (see SMS08 for details). 



4.1. Efficiency of Heat Transport 

Terrestrial planets with different rotation rates will redistribute heat from the substellar 
point (or, in a ID model, the substellar latitude) with different efficiencies. According to the 
idealized scaling described above, wherein the effective diffusion coefficient varies with the 
inverse square of the planetary rotation rate, slower spinning planets will redistribute heat 
more efficiently, while faster spinning planets will do so less efficiently. But from where, and 
to where, is heat redistributed? How does this depend on obliquity and rotation rate? And 
what influence does this have on climatic habitability? 
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4.1.1. Direction of Heat Flow 

For an Earth-like 23.5° obliquity, the substellar latitude does not vary very much over the 
course of the year: the tropics are fairly close to the equator (the tropical region is less than 
one third of the Earth's surface area). As a result, it is a reasonable approximation that 
heat is always being transported from the equator to the poles (but see Fig. [2] for details). 
In contrast, on a planet with significantly larger obliquity, the direction of heat fiow changes 
over the course of the annual cycle. At the equinoxes, the equator is the most strongly 
insolated part of the planet (regardless of obliquity), and so heat builds up at the equator, 
to be partially redistributed by atmospheric motions. But on a highly oblique planet, polar 
summers are extremely intense, as measured by diurnally averaged insolation. As a result, 
heat builds up at the poles during their corresponding summers, and the fiow of heat reverses 
direction. 

Figure [5] demonstrates the effect of such strong polar summers on the global radiation budget 
of a model planet, by comparison to Figure [T] in which the annually averaged cooling and 
heating are shown for an Earth-like, 23.5° obliquity model. As expected for the Earth-like 
model, over the annual cycle, the equator receives significantly more solar radiation than 
do the poles, and accordingly the annually averaged heating exceeds the cooling at the 
equator. This indicates that atmospheric motions transport heat poleward from the equator 
on average. In Figure [5| on the other hand, we present the analogous functions in the case of 
high and extreme obliquity models. The left panel shows the heating and cooling functions 
for a model at 60° obliquity; the right panel shows the same functions for a model at 90° 
obliquity. Planetary scientists have long recognized that in highly oblique models such as 
these, the polar summers are so intense that, averaged over the year, the most strongly 
insolated parts of the planet are the North and South Poles! (See, e.g., ?.) In an annually 
averaged sense, then, heat fiows from the poles to the equator, although clearly the direction 
of fiow changes with the seasons, as described above. The import of these plots is that 
our notion that the poles are the coldest planetary regions might have to be revised in the 
case of highly oblique worlds. The resulting regime of atmospheric transport, which is only 
parameterized with our diffusive treatment, may also be expected to differ substantially from 
that on Earth (e.g., in terms of equatorial Hadley cells). 

Figure |6] shows in greater detail the extreme way in which insolation can vary over the 
annual cycle in a highly oblique model. In this model, the obliquity of an otherwise Earth- 
like planet {I2-A2 functions, with 70% ocean uniformly distributed) is set to 90°. Notice 
that the cooling remains much more steady than the heating in this model. This is because 
of the high effective heat capacity of the atmosphere above ocean. In models with less ocean 
coverage, or oceans that are nonuniformly distributed, the cooling too can vary dramatically 
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over the annual cycle Because the heating varies so intensely while the cooling varies less, 
their difference - the net heating curve - also exhibits large variations within the annual 
cycle. At each solstice, the pole facing the star receives far more net radiant flux than 
both the opposite pole and the planetary mean. At the equinoxes, something perhaps more 
surprising happens: the net radiant flux is negative over much of the model planet's surface, 
and only barely positive near the equator. Overall, the planet heats strongly at the poles 
during solstices (while cooling elsewhere) and either cools or remains essentially thermally 
neutral everywhere during equinoxes. 



4- 1.2. Implications for Habitability 



As SMS08 demonstrates, model planets with eflicient heat transport (slowly spinning plan- 
ets) are more latitudinally isothermal than models with Earth-like rotation, which themselves 
exhibit less latitudinal variation of temperature than those with ineflicient heat transport 
(fast spinning planets). As a result, models corresponding to worlds that are spinning slowly 
relative to the Earth (but still fast enough that a 1-D climate model has some use) tend to be 
either entirely habitable or entirely non-habitable at any given time. In contrast. Earth-like 
and faster spinning model planets may be only partially habitable at a particular time. They 
may, for instance, be frozen at the poles and temperate at the equator, or vice versa in the 
case of a highly oblique world. 

Figure [7] demonstrates the complicated interplay that can go on between obliquity and ef- 
ficiency of heat transport, in determining a planet's habitability. This figure shows the 
temporal habitability fraction, as a function of orbital semimajor axis and latitude, for each 
of 12 diflFerent combinations of obliquity (0°, 30°, 60°, 90°) and latitudinal heat diflFusion 
coeflicient (i^fid/9, ^fid, and 9D^^^ corresponding respectively to 72-hour, 24-hour, 8-hour 
rotations). The top panel shows these plots for the (/2, A2) cooling-albedo pair, and the bot- 
tom panel depicts the (73,743) pair (see Table [l] for details). In both panels, the left column 
of plots represents eflicient latitudinal heat transport; the middle column represents Earth- 
like transport; and the right column represents ineflicient transport. Each of the 24 plots 
in this flgure shows the results of model runs for planets located from 0.45 AU to 1.25 AU 
in increments of 0.025 AU. The color scale indicates the fraction of the year that the lati- 
tude at that point spends in the habitable temperature range (273 K - 373 K) on a model 
planet at the specifled orbital semimajor axis. In each plot, the white vertical dashed lines 
indicate the radiative equilibrium habitable zone, calculated (as discussed in SMS08) from 



^It is also worth noting that at the outer reaches of a ~ 2Mq star's habitable zone (e.g., Kasting et al. 
1993), the annual cycle might be long enough relative to the thermal timescale (~ a decade) of the ocean- 
atmosphere mixed layer that it could undergo relatively large swings in temperature within a single annum. 
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a 0-dimensional model with annually averaged, globally averaged insolation and cooling. 

There are a number of interesting features in Figure [7| The most obvious one is that, as 
expected, at every obliquity, less efficient transport results in more strongly latitudinally dif- 
ferentiated temporal habitability. In addition, at each transport-efficiency value, the > sign 
shape of the seasonally habitable ribbon at low obliquity reverses to a < sign shape at high 
obliquity. In other words, at low obliquity, the relatively cold poles are habitable closer to 
the star and the relatively warm equator is habitable farther from the star. At high obliquity, 
however, this reverses and the poles are relatively warm, while the equator is comparatively 
colder. 



Furthermore, in both panels, the contours in most cells show a very abrupt outer boundary 
to the seasonally habitable zone. This is because, as discussed in SMS08, the ice-albedo 
feedback renders these models quite sensitive to changes in forcing. Small reductions in 
insolation can be amplified, because the ice-coverage increases, which increases the global 
albedo and leads to further reduction in insolation. This feedback mechanism renders the 
model climates susceptible to a global snowball transition, from which they cannot recover 
within our model framework. |^ The main exceptions to this trend are the low obliquity, 
fast-spinning models, in the upper right corners of both panels, although even these models 
drop to 0% habitability at orbital radii that are small relative to the outer boundary of 
the habitable zone set by global radiative equilibrium (indicated by the white dashed lines). 
Interestingly, the fairly small difference in cooling-albedo functions from (/2, A2) to (/3, ^43) is 
sufficient to allow the cell in the lower right corner of the bottom panel - extreme obliquity, 
inefficient transport - to avoid transitioning globally to a snowball state. In that model, 
the intense summer insolation at the poles, combined with the relative thermal isolation of 
different latitudes, allows the poles to heat up above the freezing point of water during their 
summers, even at orbital distances where other models would be entirely frozen. In sum, 
susceptibility to snowball transitions depends on details of parameterizations in our energy 
balance model, as has been noted before for Earth climate studies (e.g.. North et al. 1981). 



4.2. Land/Ocean Distribution 



As described in SMS08, the large covering fraction of oceans on the Earth (roughly 70%) 
stabilizes our climate over an annual cycle, by virtue of the large effective heat capacity of 
atmosphere over ocean. Over land, the thermal relaxation timescale is several months, while 



^Additional feedback mechanisms that are not incorporated in our model exist on a real planet and might 
help it to recover from a snowball state. For discussions of how the Earth might have recovered from one or 
more snowball episodes, see, for example, ?, Hoffman & Schrag (2002), and ?. 



for the 50 m deep) mixed layer over the ocean, the thermal relaxation timescale is more 
than a decade. As a result, in a ID model (such as ours) that does not resolve continents 
in longitude, any latitude band with significant ocean fraction will have strongly suppressed 
annual temperature variations relative to a latitude band with low ocean fraction. Because 
we do not know of any way to determine a priori the distribution of continents and oceans 
on an extrasolar planet, it is important to consider the infiuence on climatic habit ability of 
other possible land/ocean distributions. 



4.2.1. Nonuniform Ocean Coverage 

We consider model planets with distributions of land and ocean that are not uniform across 
different latitudes: one with 30% land coverage, with a land-mass centered on the North 



Pole (extending down to ^ 24° North latitude), and the other (discussed in § 4.2.2) with 
90% land, again centered on the North Polej^ Because of the relatively low thermal inertia 
of atmosphere over land, parts of a model planet that are dominated by land can freeze or 
boil during the course of the year and still return to temperate conditions at other times. 
In fact, at some orbital distances, and at high obliquity, the polar regions of some models 
freeze and boil within an annual cycle. 

Figure [8] displays the tremendous swings of temperature that can occur over latitude bands 
that lack ocean, and also indicates that annually averaged calculations can miss a lot of 
information about instantaneous conditions on oblique planets. This figure contrasts the 
annually averaged temperature with the detailed temperature evolution on a model planet 
with a North Polar continent that is 30% of the total surface area, at an orbital distance of 
1 AU. This model uses (/2, A2) cooling-albedo functions and results are shown for obliquities 
23.5°, 60°, and 90°. At all three obliquities, the left column - the annually averaged temper- 
ature profile ^ provides an impoverished view of the actual climatic conditions. Looking at 
just the left panels: the 23.5° obliquity model appears slightly asymmetrical in temperature 
distribution, with the continental North Pole 8 K warmer than the oceanic South Pole; the 
60° obliquity model appears cooler at the continental pole; and the 90° obliquity model again 
appears warmer at the continental pole, but appears frozen over the whole globe. In truth, 
all three models reach significantly higher temperatures at the continental pole during its 
summer than at the other pole. At both 60° and 90° obliquity. North Pole summer tempera- 
tures exceed 410 K, as the Sun shines nearly straight down on the pole for months. We note 
that an important limitation of our models is apparent in this figure. Although we may not 
have much intuition for what the polar summers should be like on high obliquity planets, 
it is surprising to obtain summer polar continent temperatures in excess of 310 K in the 



^ Equivalent ly, this model can be conceived as having an ocean centered at the South Pole. 
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23.5° obliquity model. Indeed, Antarctica - Earth's continental pole - is significantly colder 
than the non-continental pole, and for the most part neither pole ever reaches temperatures 
above freezing. Accounting for Antarctica in our model framework would be non-trivial, as it 
may require initial conditions diflFering from uniformly ice- free and/or improved treatments 
of ice-related surface processes. 

Figure [9] presents plots of the temporal habitability fraction for the same model planet as 
above, at obliquities 0°, 30°, 60°, 90°. Figure [8] illustrated how the presence of land at the 
North Pole causes tremendous swings in temperature there; Fig. [9] confirms that, for nonzero 
obliquity, this is indeed the case throughout the orbital extent of the habitable zone. These 
large seasonal variations lead to exotic shapes in plots of temporal habitability. Compared 
with uniformly ocean-dominated worlds, much more of the parameter space is partially 
habitable at each obliquity (except 0°), neither 0% nor 100% of the year, but somewhere in 
between. 



4.2.2. Desert Worlds 

We now consider two model planets with just 10% ocean fraction. We examine the cases of 
a uniformly distributed ocean (10% in every latitude band) and an ocean localized at the 



South Pole (extending northward to ^ 53° South latitude). Figure 10 presents the temporal 
habitability for the uniform desert world, while Figure [TT] presents the analogous plot for 
the desert world with a South Polar ocean. As before, some regions of these model planets 
swing from freezing to boiling temperatures over the course of the year. This is responsible 
for the butterfly shape of the temporal habitability plots in the 60° and 90° obliquity cases 



shown in Fig. 10: at a 0.9 AU, the poles are habitable for a smaller fraction of the year 



than more equatorial regions at that orbital distance, or than the poles for closer and more 



distant orbits. The pattern of habitability in Fig. 11, on the other hand, shows how strongly 



assymmetric the climate can be on a desert world with a polar ocean. 



These models, and those presented in § |4.2.1[ suggest that at extreme obliquity the inner 



edge of the zone of regionally and seasonally habitable climates can be extended dramatically 
inward, while the outer boundary can only be extended mildly outward. Several important 
caveats should accompany this observation, however. Assuming an infrared cooling function 
that is constant with orbital radius probably leads to a flawed treatment at both the high and 
low insolation limits of these models. At the inner edge of the habitable zone, large increases 
in atmospheric water content can cause a reduction in the cooling efliciency, leading to a 
runaway greenhouse eflFect. An eventual catastrophic water loss can result in a Venus-like 



outcome, as described by Kasting et al. (1993) and references therein (although this type of 



outcome might be mitigated by reduced heating from increased cloud-albedo, as mentioned 
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in SMS08). At the outer edge, on long timescales, a reduced efficiency of the carbonate- 
sihcate weathering cycle is likely to lead to a significant increase in the partial pressure of 
atmospheric CO2 ( Kasting et al.|p!993 ), which could extend the habitable zone from ^ 1 AU, 



as in our models, to ^ 1.4 AU or beyond in some cases. We discuss in greater detail the 



issue of varying atmospheric CO2 content with orbital distance in § 4.3 



Notwithstanding these various complications, for plausible cooling functions, the low ther- 
mal inertia of atmosphere over land might lead to severe polar climates on highly oblique 
planets. What are we to make of partial "habit ability" by our criterion in the case of a 
region of a planet that actually boils and freezes every year? There are some microbes on 
Earth that can reproduce at freezing temperatures, and others that can reproduce at boiling 
temperatures (see SMS08 and references therein) , although none of which we are aware that 
can do both. If part of a planet regularly swings through these wild extremes of climate, 
is it appropriate to call it habitable? More relevant from the perspective of formulating a 
testable scientific hypothesis, could such a planet support enough life to produce sufficient 
levels of biosignatures that its life could be detected from Earth? This is an open question, 
but it is worthwhile to keep in mind that microbes on Earth appear to be as hardy as they 
need to be: nearly everywhere that biologists have searched, they have found some microbes 
thriving. A perhaps equally significant result from the perspective of habitability is that the 
reduced thermal inertia of these models appears to render them somewhat less susceptible 
to global snowball events, especially at high obliquities (e.g., compare Figs. [9 11 with the 
middle column of Fig. [7|. 



4.3. Modeling the Far Reaches of the Habitable Zone 



Walker et al. (1981D propose that a planet's temperature is regulated on long timescales by a 



feedback mechanism involving weathering of silicate rocks through carbonic acid from CO2 
disolved in water. They argue that since the rate of weathering (and hence of removal of 
CO2 from the atmosphere) increases with temperature, this process is an important negative 
feedback on climate that acts to keep temperatures near the freezing point of water (see the 
recent study by ? confirming the operation of this cycle). Kasting et al. ( 1993[ ) point 
out that this negative feedback can significantly offset the extreme sensitivity of climate to 



changes in orbital distance away from 1 AU seen in models such as those of Hart (1979) 



and the models presented in SMS08 and thus far in this paper. These models are sensitive 
because they contain a significant positive feedback of the Earth-climate - the ice-albedo 
feedback whereby at lower temperatures the absorbed insolation is dramatically reduced 
because of the high albedo of ice. However, these models also ignore the long-term, negative 
(counterbalancing) feedback of the carbonate-silicate weathering cycle ( [Kasting et al.|[T993 ). 
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In order to probe the combined influence on climate of rotation rate and obliquity in the 
context of the expected C02-rich atmosphere that a pseudo Earth would have at 1.4 AU, we 
switch to the infrared cooling function used by WK97 (/wk97 in Table [T]), with CO2 partial 
pressure {PCO2) set to 1 bar and 2 bars. For simplicity, we maintain the same albedo function 
A2 as before and adopt the same linear dependence of the latitudinal heat diffusion coeflicient 
with total atmospheric pressure as WK97 {D/D^^ oc Ptot)- We flnd that at both 1- and 2- 
bar levels of CO2, model planets maintain globally temperate conditions at all obliquities for 
both D = 9Dfid and D = ^ corresponding to slow and Earth-like rotation. Interestingly, 
however, at reduced transport efliciency, corresponding to fast planetary rotation, we flnd 
that these C02-rich models are susceptible to global glaciation events. 



Figure [12] shows the global average temperature and the climate evolution for fast-spinning 
model planets at 1.4 AU with 1 bar atmospheric CO2, (/wk97,^2) cooling-albedo functions 
and a latitudinal heat diffusion coefflcient scaled up with total pressure but reduced by a 
factor 9 to account for rapid rotation. At 23.5° obliquity, the cold temperatures at the poles 
drag the model into a snowball state, while at 90° obliquity, it is the cold equator that 
drags the model planet to the same fate. At 60° obliquity, however, no part of the planet 
receives consistently low enough insolation to trigger global glaciation. The snowball effect 
seen in the 23.5° and 90° obliquity models might be particularly dramatic because of the 
possibility of partial atmospheric collapse, on a much shorter timescale than volcanism can 
replenish CO2. At 1 bar, the freezing point of CO2 is - 195 K. Both the 23.5° and 90° 
obliquity models reach temperatures below this threshold over large enough regions of their 
surface^ that signiflcant amounts of the atmospheric CO2 might condense out as dry ice, 
thereby reducing the atmospheric greenhouse effect. The risk of atmospheric collapse would 
be somewhat lessened by the release of latent heat from CO2 condensation, which would tend 
to prevent too much CO2 from freezing out during any winter. Partial collapse, on the other 
hand, would reduce surface pressure and thus the efflciency of atmospheric heat transport. 
Realistically treating these possibilities is beyond the scope of the present paper as it would 
require incorporating a latent heat term in the energy balance equation and accounting for 



the surface CO2 mass budget, as Nakamura & Tajika (2002) do in their Mars EBM. 



While model runs with 2 bars of CO2 (not shown) indicate that a pseudo-Earth with such a 
massive atmosphere at 1.4 AU would be unlikely to suffer glaciation or partial atmospheric 
collapse at any obliquity, it is worth noting that on a planet with a rate of volcanic green- 
house gas output not much higher than on Earth, it may be difflcult to build up a thick CO2 
atmosphere (e.g., 2 bars) under conditions such that CO2 condenses out at lower pressure 
(like in our model with 1 bar of CO2). In situations in which atmospheric CO2 concentrations 
can change appreciably on a yearly timescale, our model is not self-consistent. Further anal- 



^The very large seasonal variations of temperature on these snowball planets result from the moderate 
heat capacity of the atmosphere over ice (see, e.g., SMS08; WK97). 
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ysis is therefore needed to determine the fuU extent and consequences of such atmospheric 
condensation events. Clearly, these various issues are important subjects for future stud- 
ies since they indicate that snowball transitions could in principle limit the habitability of 
some terrestrial planets by interfering with the negative feedback of their carbonate-silicate 
weathering cycles. 



4.4. Validity of Global Radiative Balance 



Calculations of habitable zones have often assumed global radiative balance conditions. Al- 
though these calculations by definition cannot account for the regional character of habitabil- 
ity, one might hope that they still provide a decent proxy for the global average conditions. 
Indeed, as seen in Fig. |3} the Earth itself is within ^ 5% of radiative equilibrium throughout 
the year. Accordingly, global radiative balance models have provided a very useful starting 
point for considerations of how the habitability of an Earth-like planet depends on its orbital 
radius. 



Figure [13] presents globally averaged cooling, heating, and net radiative fiuxes in a model 
with a North Polar continent that covers 30% of the planet's surface, with 90° obliquity. 
In contrast to the Earth, which remains within 5% of global radiative balance, this model 
planet can be nearly 60% out of global radiative balance. While it has been recognized 
that planets on highly eccentric orbits experience forcings that are significantly different 
from annually and globally averaged conditions, the results in Fig. [13] illustrate how, even 
on circular orbits, planets can experience conditions that are far from radiative equilibrium. 
This further underscores the importance of regional, time-dependent climate models for 
addressing the habitability of extrasolar terrestrial planets. 



5. Conclusions 



We have presented a series of energy balance models to address the variety of climatic con- 
ditions that might exist on oblique terrestrial planets with circular orbits. We considered 
dynamic climate forcings and responses determined by several planetary attributes a priori 
unknown for extrasolar planets, including obliquity, rotation rate, distribution of land/ocean 
coverage, and the detailed nature of the radiative cooling and heating functions. We find that 
planets with small ocean fractions or polar continents can experience very severe seasonal 
climatic variations, but that these planets also might maintain seasonally and regionally hab- 
itable conditions over a larger range of orbital radii than more Earth-like planets. Climates 
on high obliquity planets with nonuniform distributions of land and ocean can be far from 
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global radiative balance, as compared to the Earth. Our results provide indications that the 
modeled climates are somewhat less prone to dynamical snowball transitions at high obliq- 
uity. Fast rotating Earth-like planets may fall victim to global glaciation events at closer 
orbital radii than slower rotating planets. This is also the case for planets with massive CO2 
atmospheres, which are expected to be found in the outer orbital range of habitable zones. 
Snowball transitions could be particularly significant for such planets since partial collapse 
of their C02-rich atmospheres may occur and possibly interfere with the thermostatic effect 
of their carbonate-silicate weathering cycle, thus affecting their long-term habitability. 
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Table 1: Atmospheric Models. 



Model IR Cooling Function Albedo Function 

h[T] = i+(3;;;a[T T ^2[r]= 0.525 - 0.245 tanh[i;^^] 
3^ h [T]=A + BT [T] = 0.475 - 0.225 tanh[^^^^] 

4^ JwK97[r,pC02] ^[r] 



Model with T-dependent optical thickness: TiR[r] = 0.79(r/273K)3. 
''Linearized model: A = 2.033 x 10^ erg cm~^ s~\ B = 2.094 x 10^ erg cm~^ s~^ 

"^WKQ? cooling function (with A2 albedo). The detailed functional form is presented in the appendix of 
WK97. 

Note. — a is the Stefan-Boltzmann constant. 
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Fig. 1. — Annually averaged cooling and heating fluxes for our flducial model at 1 AU and 
for the Earth. The thick red line is infrared cooling and the thick blue line is absorbed solar 
flux, in our flducial model (70% uniform ocean, /2, A2). The thin dashed magenta line is the 
Earth's annually averaged long wavelength infrared radiation, and the thin dashed cyan line 
is the annually averaged absorbed solar flux measured on Earth. Earth-speciflc features are 
not captured by our symmetric and uniform model. 
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Fig. 2. — Monthly cooling, heating, and net fluxes for our flducial model at 23.5° obliquity 
and 1 AU orbital distance, and for the Earth. Each panel presents the average cooling, heat- 
ing, and net (heating minus cooling) radiative fluxes, as functions of latitude, for one month 
of the year, starting at the Northern winter solstice (upper left panel), and incrementing by 
one month with each panel to the right. These fluxes are presented for both the model (thick 
solid lines) and the Earth (thin dashed lines). Model heating is blue; model cooling is red; 
model net heating is green. Earth heating is cyan; Earth cooling is deep magenta; Earth net 
heating is black. Our model captures reasonably well the seasonal variations of these fluxes. 
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Fig. 3. — Global average cooling, heating, and net radiative flux, as functions of time in 
our flducial model at 23.5° obliquity (solid lines). Earth's ERBE data are shown as stars. 
Top Panel: Global average cooling (red curves and magenta stars) and heating (blue curves 
and cyan stars) fluxes as a function of time of year, measured in fraction of a year from the 
Northern winter solstice. Bottom Panel: Net heating flux (heating minus cooling) for the 
model (green curve) and ERBE data (black stars), plotted as percent of the corresponding 
cooling flux. The Earth remains within 5% of global radiative balance throughout its seasonal 
cycle. 
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Fig. 4. — Annually averaged, longitudinally integrated, meridional heat transport rates in 
models at 23.5° and 90° obliquity. Transport is positive northward. In the Earth-like case 
(blue solid curve), heat flows from the equator to the poles; in the highly oblique case (green 
dashed curve), the annually averaged heat flow is reduced and in the opposite direction. This 
result is in close agreement with a comparable one obtained by Williams & Pollard (2003) 
with a full physics climate model. 




Fig. 5. — Annually averaged cooling and heating fluxes for high and extreme obliquity model 
planets, compared to Earth. By contrast with Fig. [T} which shows the annually averaged 
heating and cooling for 23.5° obliquity, here we present these functions for Earth-like models 
at 60° obliquity (left panel) and at 90° obliquity (right panel). Thick lines are model results; 
thin dashed lines are Earth's ERBE data. In both high obliquity cases, unlike on Earth, there 
is net annually averaged heating at the poles (i.e., heating exceeds cooling), and, especially 
for the extreme obliquity case, net annually averaged cooling at the equator. 




Fig. 6. — Monthly cooling, heating, and net heating fluxes for the flducial model at 90° 
obliquity and 1 AU orbital distance, and for the Earth. Each panel presents the average 
cooling, heating, and net (heating minus cooling) radiative fluxes, as functions of latitude, 
for one month of the year, starting at the Northern winter solstice (upper left panel), and 
incrementing by one month with each panel to the right. These fluxes are shown for both the 
oblique model planet (thick solid lines) and the Earth (thin dashed lines). Model heating is 
blue; model cooling is red; model net heating is green. Earth heating is cyan; Earth cooling 
is deep magenta; Earth net heating is black. While cooling fluxes remain relatively steady 
and uniform on such an oblique planet, heating and net heating fluxes are subject to very 
large seasonal variations. 
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Fig. 7. — Temporal habitability fraction in models with different obliquities, heat transport efficiencies, 
and cooling-albedo functions. In both panels, obliquity varies from 0° (top row) to 90° (bottom row). The 
latitudinal heat diffusion coefficient varies from (left column) to I^fid (center) and 9I)fid (right). In 

each panel, the abscissa is orbital radius, in AU, and the ordinate is latitude. Colors indicate the fraction 
of the year spent by that region in the habitable temperature range (273 K - 373 K). Top Panel: (12,^2) 
cooling-albedo combination. Bottom Panel: (73,743) cooling-albedo combination. For comparison, vertical 
dashed lines show the habitable zone extent expected on the basis of global radiative balance. Regionality 
and seasonality can extend the inner reach of the instantaneous habitable zone while often reducing its outer 
reach, when global snowball transitions occur. 
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Fig. 8. — Annually averaged and detailed time-dependent temperature profiles in models with a North 
Polar continent that takes up 30% of surface area (the other 70% is ocean) at 1 AU orbital distance, for 
(12,^2) cooling-albedo functions and obliquities of 23.5°, 60° and 90°. Left: The magenta curve shows the 
annually averaged temperature profile for models with a North Polar continent that extends down to ~ 24° 
North latitude. The solid blue and dashed green curves are for reference - blue: the fiducial model (identical 
to this one, except for 70% ocean uniformly distributed in every latitude band); dashed green: the Earth's 
actual temperature profile, as measured by NCEP/NCAR in 2004 ( [Kistler et ar]|1999[ [Kalnay et al.||1996 ). 
Right: Detailed temperature evolution (latitude vs. time) in the polar continent models from years 8 through 
50. Top Row: 23.5°. Middle Row: 60°. Bottom Row: 90°. Annually averaged profiles (left) miss much 
of the seasonal variations shown by the detailed profiles (right). The 90° obliquity model experiences an 
asymmetric, partial snowball transition. 
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Fig. 9. — Temporal habitability fraction at different obliquities for a model planet with a 
North Polar continent covering 30% of its surface. Obliquity varies from 0° to 90° from 
top to bottom. The latitudinal heat diffusion coefficient is kept to its fiducial value. The 
notation is similar to Figure [7} colors indicate the fraction of the year spent by that region 
in the habitable temperature range (273 K - 373 K). Habitability has a strong regional and 
seasonal character on such a planet, when oblique. 
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Fig. 10. — Temporal habitability fraction at different obliquities for a model planet with 
10% ocean uniformly distributed. Otherwise similar to Figure |9j Habitability has a strong 
regional and seasonal character on such a planet, when oblique. 
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Fig. 11. — Temporal habitability fraction at different obliquities for a model planet with 
a North Polar continent covering 90% of its surface (i.e., a localized South Polar ocean). 
Otherwise similar to Figure [9| Habitability has a strong regional and seasonal character on 
such a planet, when oblique. 
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Fig. 12. — Annually averaged and detailed time-dependent temperatures on a fast spinning 
model planet at 1.4 AU, using the infrared cooling function of WK97 for a C02-ricli atmo- 
sphere with J9CO2 = 1 bar. From top to bottom, results are shown for an obliquity of 23.5°, 
60° and 90°. The notation and the models are the same as in Figure [8} except for the cooling 
function (/wk97 with PCO2 = 1 bar) and a latitudinal heat diffusion coefficient adjusted for 
fast spin and a massive atmosphere. In the two globally- frozen models (top and bottom), 
surface temperatures that would result in seasonal CO2 atmospheric collapse (T ^ 195 K) 
are reached over a fair fraction of the planet's surface area. 
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Fig. 13. — Global average cooling, heating, and net (heating — cooling) radiative flux, as 
functions of time, for a model planet with a North Polar continent covering 30% of its surface, 
at 90° obliquity. Notation is similar to Figure [3} This planet experiences large deviations 
from global radiative balance. 
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